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AN ELEMENTARY SOLUTION OF GESSEL’S WALKS 
IN THE QUADRANT 

MIREILLE BOUSQUET-MELOU 


Abstract. Around 2000, Ira Gessel conjectured that the number of lattice 
walks in the quadrant N 2 , starting and ending at the origin (0,0) and taking 
their steps in {—>, /A •<—, had a simple hypergeometric form. In the follow¬ 
ing decade, this problem was recast in the systematic study of walks with small 
steps (that is, steps in {—1,0, l} 2 ) confined to the quadrant. The generating 
functions of such walks are archetypal solutions of partial discrete differential 
equations. 

A complete classification of quadrant walks according to the nature of their 
generating function (algebraic, D-finite or not) is now available, but Gessel’s 
walks remained mysterious because they were the only model among the 23 
D-finite ones that had not been given an elementary solution. Instead, Gessel’s 
conjecture was first proved using an inventive computer algebra approach in 
2008. A year later, the associated three-variate generating function was proved 
to be algebraic by a computer algebra tour de force. This was re-proved re¬ 
cently using elaborate complex analysis machinery. We give here an elementary 
and constructive proof. Our approach also solves other quadrant models (with 
multiple steps) recently proved to be algebraic via computer algebra. 


1. Introduction 

The enumeration of planar lattice walks confined to the quadrant has received 
a lot of attention in the past decade. The basic question reads as follows: given 
a finite step set S C Z 2 and a starting point F e N 2 , what is the number q(n) 
of n-step walks, starting from P and taking their steps in S , that remain in the 
non-negative quadrant N 2 ? This is a generic and versatile question, since such 
walks encode in a natural fashion many discrete objects (systems of queues, Young 
tableaux and permutations among others). More generally, the study of these walks 
fits in the larger framework of walks confined to cones. They are also much studied 
in probability theory, both in a discrete [22, 23] and a continuous [21, 27] setting. 
From a technical point of view, counting walks in the quadrant is part of a general 
program aiming at solving partial discrete differential equations , see (5) and (7) 
below. Details are given in Section 2, together with a more algebraic viewpoint 
involving division of power series. 

On the combinatorics side, much attention has focused on the nature of the 
associated generating function Q(t) = '}2 n q(n)t n . Is it rational in t, as for un¬ 
constrained walks? Is it algebraic over Q (f), as for walks confined to a (rational) 
half-space? More generally, is it D-finite, that is, a solution of a linear differential 
equation with polynomial coefficients? The answer depends on the step set, and, 
to a lesser extent, on the starting point. 

A systematic study was initiated in [16, 40] for walks starting at the origin 
(0,0) and taking only small steps (that is, S C {—1,0, l} 2 ). For these walks, a 
complete classification is now available. In particular, the generating function Q(t) 
(or rather, its three-variate refinement Q(x,y;t) that also records the coordinates 
of the endpoint of the walk) is D-finite if and only if a certain group of rational 
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transformations is finite. The proof involves an attractive variety of tools, ranging 
from basic power series algebra [12, 16, 40] to complex analysis [38, 45], computer 
algebra [7, 31] and number theory [10]. 

In the D-finite class lie, up to symmetries, exactly 23 step sets S (often called 
models ) with small steps. A uniform approach, the algebraic kernel method , es¬ 
tablishes the D-finiteness of 19 of them [16]. These 19 models are transcendental. 
The remaining 4 models are structurally simpler, since they are algebraic, but they 
are also harder to solve. Three of them, including the so-called Kreweras’ model 
S = {•<—,4-,/'}, were solved in a uniform manner in [16], following several ad hoc 
proofs [35, 28, 12, 11], The final one, S = {-», /*,«—, j/}, has been recognized as 
extremely challenging, and has been promoted recently as a popular model in the 
French magazine Pour la science [9], It is named after Ira Gessel, who conjectured, 
around 2000, that for this model the number of walks of length 2 n ending at the 
origin is 


(5/6) w (l/2) ra 
(5/3)„(2) 


where (a) n = a(a + 1) • • • (a + n — 1) is the ascending factorial. This is sequence 
A135404 is Sloane’s Encyclopedia of integer sequences [30]. The associated gen¬ 
erating function is algebraic, but no one noticed it at that time. This attractive 
conjecture remained open for several years. In particular, it resisted the otherwise 
complete solution of D-finite models presented in [16]. Finally, Kauers, Koutschan 
and Zeilberger [31] succeeded in proving (1) via a clever computer-aided approach. 
About a year later, this was generalized by Bostan and Kauers [7], who proved, by 
a tour de force in computer algebra, that the complete generating function Q(x, y\ t) 
is also algebraic. This series is defined as 


Q(x,y;t)= ^2 q{i,j;n)x l y J t n , (2) 

i,j,n>0 


where q(i,j;n) is the number of n-step walks in the quarter plane that start from 
(0,0), end at ( i,j) and take their steps in S. A complete description, and ratio¬ 
nal parametrization of this series is given in Theorem 1 below. Both the result 
of [7] and the immense computational effort it required were remarkable, but they 
left open the problem of finding a “human” proof of Gessel’s conjecture, and of its 
trivariate generalization. More recently, Bostan, Kurkova and Raschel [8] proposed 
another proof, which involves the deep machinery of complex analysis first devel¬ 
oped for stationary random walks [23] , and subsequently for counting walks in the 
quadrant [38, 45, 24], The solution involves Weierstrass’ p and £ functions, which 
are transcendental, and may appear a surprisingly complicated detour to establish 
algebraicity. Let us also cite a number of other attempts [1, 24, 33, 37, 42, 45, 47]. 
Some of them can be seen as preliminary steps to the proofs of [31] or [8]. 


In this paper, we present the first elementary solution of Gessel’s walks in the 
quadrant. It is elementary in the sense that it remains at the level of formal 
power series and polynomial equations. Also, it involves no guessing nor high- 
level computer algebra 1 . Finally, it is constructive, in the sense that we construct, 
incrementally, the extensions of Q(x, y, t) that are needed to describe the series 
Q(x,y;t). Even if this sounds a bit technical at this stage, we can inform the 
quadrant experts that the key ingredient is that the symmetric functions of the 
roots of the kernel are polynomials in 1/x. This is also the case for Kreweras’ 
walks, and our proof is in this sense close to one solution of Kreweras’ model [11]. 


1 Although one feels better with a computer at hand when it comes to handling a system of 
three polynomial equations. An accompanying Maple session is available on the author’s webpage. 
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Theorem 1 . The generating function Q(x,y,t) is algebraic over Q(x, y 1 1), of de¬ 
gree 72. The specialization Q(0,0;f) has degree 8, and can be written as 


<2(0,0; i) 


32Z 3 {3 + 3Z -3Z 2 + Z 3 ) 
(l + Z)(Z 2 + 3 ) 3 


where Z = \/T and T is the only power series in t with constant term 1 satisfying 

r= 1 + 256i2 pW (3) 

The series Q(xt, 0; t ) is an even series in t, with coefficients in Q[x], and is cubic 
over Q(Z,x). It can be written as 

16 T(U + UT — 2T)M(U, Z) 


Q(xt,0;t) = 


(1 - T)(T + 3) 3 (U + Z){U 2 -9 T+ 8 TU + T 2 - TU 2 ) 


where 


M(U, Z) = (T — 1 ) 2 U 3 + Z(T - 1 )(T - IQZ - T)U 2 

- TU(T 2 + 16 ZT - 82 T - 16Z + 17) - ZT(T 2 - 18 T + 128Z + 81), 

and U is the only power series in t with constant term 1 (and coefficients in Q[:r],) 
satisfying 

16 T 2 (U 2 -T) = x(U + UT- 2T){U 2 - 9T + 8 TU + T 2 - TU 2 ). 

Finally, the series Q(0,y;t) is cubic overQ(Z,y). It can be written as 

1&VZ 3 (3 + V+ T-VT)N(V,Z) 

^ ,J/ ’ (T-1)(T + 3) 3 (1 + V) 2 (1 + Z +V -VZ) 2 

where 

N(V, Z) = (Z- 1 ) 2 (T + 3)V 3 + (T — T)(T + 2Z - 7)V 2 

+ (T- 1)(T -2 Z- 7)V + {Z + T) 2 (T + 3), 

and V is the only series in t, with constant term 0, satisfying 

1-T + 3V + VT = yV 2 {3 + V + T- VT). 

The algebraicity was established by Bostan and Kauers, and the parametrization 
was given by van Hoeij in the appendix of their paper [7] (with his notation, T = v, 
U = u and V = —w). The parametrization that we construct through our proof is 
actually a bit different (and of course equivalent), but we prefer to give van Hoeij’s 
to avoid confusion. Just for the record, ours reads as follows: 

T = t 2 (l — T)(l + 3T) 3 , Z = f(l- Z + Z 2 ), (4) 

~ r 1 ( ~~(U + Z -UZ + UZ 2 )(U -UZ ~ Z 3 + UZ 2 ) 

+ X U(Z-l)(Z + l) 3 (Z 2 + U) ’ 

- {V-1) 2 {1-Z-V + Z 2 ~ Z 3 ) 

V = y -=-=-=-, 

(1 -Z + V-Z 2 ) 2 

with constant terms 0 for T and Z, 1 for U and y/( 1 + y) for V. 

We conclude this introduction with some notation. For a ring R , we denote by 
i?[ir] the ring of polynomials in x with coefficients in R. If R is a field, then R(x) 
stands for the field of rational functions in x. Finally, if F(x\t) is a power series 
in t whose coefficients are Laurent polynomials in x, say 

F{x]t)= ^2 /(*; n)t n x l , 

n>0, 
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we denote by \x-]F(x\t) the non-negative part of F in x: 

[x-]F(x]t):= Y f(i\n)t n x\ 

n> 0, iGN 

This generalizes the standard notation \x l ]F(x',t) for the coefficient of x 1 in F. We 
similarly define the non-positive part of F(x m ,t), denoted by [x-]F(x;1). 


2. Context: partial discrete differential equations 


The starting point of the systematic approach to the enumeration of quadrant 
walks is a functional equation that characterizes their complete generating function 
Q(x,y;t ) = Q(x,y ), defined as in (2). For instance, for square lattice walks (with 
steps f,«—,4-) starting at (0,0), this equation reads: 


Q(x,y) 


1 + t(x + y)Q(x,y) + t 


Q{x,y) - Q(0,y) 

X 


t Q(x,y) - Q(x, 0 ) 

y 


( 5 ) 


It translates the fact that, to construct a walk of length n > 1, we simply add a step 
to a shorter walk of length n — 1. The divided differences (or discrete derivatives) 


Q(x,y) - Q(0,y) 

X 


and 


Q{x,y) - Q(z,Q) 
y 


( 6 ) 


arise from the fact that one cannot add a West step if the shorter walk ends on 
the y-axis, nor a South step if it ends on the x-axis. For Gessel’s walks, with steps 
-A, /*,•<—, yf' , the corresponding equation reads 


Q(x, y) = l + t(x + xy)Q(x , y)+t 1 —— —- 

X 

[ t Q(x, y) - Q(x, 0) - Q( 0, y) + Q(0,0) ^ 

xy 

because the ^ step cannot be appended to a walk ending on the x- or y-axis. Refer- 
ing to terms like (6) as discrete derivatives, it makes sense to call such equations 
partial discrete differential equations. 

The more algebraically inclined reader will write (5) as 

xy = (xy - t(x 2 y + xy 2 + x + y)) Q(x, y) + txQ(x, 0) + tyQ( 0, y), 


and observe that we are just performing the formal division of xy by the polyno¬ 
mial (xy — t(x 2 y + xy 2 + x + y)) with initial monomial xy (this is an instance of 
Grauert-Hironaka-Galligo division, see e.g. [46, Thm. 10.1]). 

Discrete differential equations are ubiquitous in combinatorial enumeration. For 
instance, for the much simpler problem of counting walks with steps ±1 starting 
at 0 and remaining at a non-negative level, the corresponding functional equation 
reads 

F(x) = F(t r; t) = 1 + txF(x) + t ~ F ^ , 

x 

an ordinary discrete differential equation (since the derivatives are only taken with 
respect to x). More generally, the enumeration of lattice walks confined to a ratio¬ 
nal half-space systematically yields linear ordinary discrete differential equations 
(DDEs) [3, 17]. Note that the order of the equation increases with the size of the 
down steps: if, in the above problem, we replace the step —1 by a step —2, the 
equation reads 


F(x) 


1 + txF(x) + t 


F(x)-F(0)-xF’(0) 

x 2 
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and involves a discrete derivative of order 2. Moving beyond linear equations, 
the enumeration of maps (connected graphs properly embedded in a surface of pre¬ 
scribed genus) has produced since the sixties a rich collection of non-linear ordinary 
discrete differential equations (see e.g. [4, 20, 48]). For instance, the enumeration 
of planar maps is governed by the following equation: 

F(z; t) = F(x) = 1 + tx 2 F(x) 2 + tx xF ( x ) ~ _ 

x — 1 

(Note that the discrete derivative is now taken at x = 1.) 

An approach for solving linear ordinary DDEs was initiated by Knuth in the 
seventies [34, Section 2.2.1, Ex. 4] and is now known as the kernel method [2, 17]. 
Another approach, developed by Brown in the sixties and known as the quadratic 
method , solves quadratic ordinary DDEs of first order [19, 29]. More recently, this 
approach was generalized by the author and A. Jehanne [15] to ordinary DDEs of 
any degree and order, with the striking result that any series solution of a (well- 
founded) equation of this type is algebraic 2 . This result is of particular importance 
to this paper, since the core of our proof is to derive from the partial DDE (7) 
satisfied by Q(x,y) an ordinary DDE satisfied by Q{x, 0). This is achieved in 
Section 3.3. This ordinary DDE (of order 3 and degree 3) is then solved using the 
method of [15], and found, of course, to have an algebraic solution. 

If ordinary DDEs are now well-understood, much less is known about partial 
DDEs, even linear ones. The enumeration of quadrant walks with small steps 
appears as an essential step in their study, because the resulting equations are as 
simple as possible: linear, first order, and involving derivatives with respect to 
two variables only. This study started around 2000, and it was soon understood 
that the solutions would not always be algebraic, nor even D-finite [18, 41]. The 
classification of these problems is now complete. That is, one knows for every 
step set S C {—1,0, l} 2 if the series Q(x,y;t) is rational, algebraic, D-finite or 
not [7, 10, 16, 38, 41, 39]. A crucial role is played by a certain group of rational 
transformations associated with S: the series Q(x,y;t) turns out to be D-finite if 
and only if this group is finite. In this case, a generic method solves the partial 
DDE associated with the problem [16]. except for four sets S for which the series 
Q(x, y ; t) is in fact algebraic. Three of these four problems were solved in a uniform 
way in [16]. The fourth one is Gessel’s model, the story of which we have told 
in the introduction. This paper thus adds a building block to the classification of 
quadrant walks, by deriving from the partial DDE (7), in a constructive manner, 
an algebraic equation for Q(x,y). 

To finish, let us mention a few investigations on more general partial DDEs: 
quadrant walks with larger backwards steps yield linear equations of larger or¬ 
der [25]; coloured planar maps yield non-linear equations [5, 14]; walks confined to 
the non-negative octant N 3 yield linear equations with three differentiation vari¬ 
ables [6]; the famous enumeration of permutations with bounded ascending sub- 
squences (and in fact, many problems related to Young tableaux, plane partitions 
or alternating sign matrices [49]) yields linear equations with an arbitrary number 
of differentiation variables [13]. 


3. The proof 

We begin, in a standard way, by writing a functional equation satisfied by 
Q(x , y\ t). It is based on a step by step construction of walks [16]: 

Q(x , y) = 1 + t(x + xy + x + xy)Q(x , y) -tx( 1 + y)Q( 0, y) - txy (Q(x, 0) - Q(0,0)), 

2 As learnt recently by the author, this algebraicity also follows from a difficult theorem on 

Artin’s approximation with nested conditions , see [36] or [43, Thm. 1.4]. 
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where Q(x,y) = Q(x,y;t), x = 1/x and y = 1/y. Equivalently, 

xyK(x, y)Q{x , y) = xy - t{Q( x, 0 ) - Q( 0 , 0 )) - t{ 1 + y)Q{ 0 , y), (8) 

where 

K (x, y) = 1 — f (x + xy + x + xy) 

is the kernel of the equation. This Laurent polynomial is left invariant by the 
following two rational transformations: 

$ : (x, y) (xy, y) and T := (x, y) (x, x 2 y). 

Both are involutions, and they generate a group G of order 8: 

(x, y)<A-(xy, y)*^->{xy, x 2 y)<A>(x, x 2 y)e-^(x, y)<A- 

{xy, y)^(xy, x 2 y)^(x, x 2 y)^(x, y). 

The construction of this group is also standard (see for instance [16]). 


3.1. Canceling the kernel 

As a polynomial in y, the kernel A'(x, y) has two roots, which are series in t with 
coefficients in Q[x,x]: 


Y 0 (x) 


1 — t(x + x) — \/{l — t(x + x)) 2 — 4f 2 
2 tx 


xt + 0(t 2 ), 


Yi{x) 


1 — t{x + x) + y/(l — t(x + x)) 2 — 4t 2 
2 tx 


|-(1+ x 2 ) — xt + 0{t 2 ). 


Observe that the series xY-i(x) are symmetric in x and x: 


xYi(x) = xYi{x). 


Moreover, the elementary symmetric functions of the Yj. namely 

Yo + Y-l = — 1 + ^ — x 2 and Fo^i = S 2 , (9) 

are polynomials in x = 1/x. This property also holds for Kreweras walks [11, 12], 
and plays a crucial role in our proof. The following lemma tells us how to extract 
the constant term of a symmetric polynomial in Yq and Y\. 


Lemma 2. Let P(u,v ) be a symmetric polynomial in u and v with coefficients 
in Q. Then P{Yq,Y\) is a polynomial in x with coefficients in Q[l/t]. Its constant 
term in x (equivalently, its non-negative part in x) is P( 0, —1). 

Proof. The first statement follows from (9) and the fact that every symmetric poly¬ 
nomial in u and v is a polynomial in u+v and uv. By linearity, it suffices to check the 
second statement when P(u,v) = u m v n + u n v m . If min(m,n) > 0, then P{Yq,Yi) 
has a factor YqY] = x 2 and its constant term is 0. If P(u,v) = u n + v n , then 
one proves by induction on n that the constant term is 2 if n = 0, and (—1)™ 
otherwise. □ 


Application. We will typically apply this lemma to prove that if F(u,w,t) is a 
series in t whose coefficients are symmetric polynomials in u and v, with rational 
coefficients, and if F{Yq, Y\ ; f) is well-defined as a series in t (with coefficients in 
Q[x,x]), then its coefficients actually lie in Q[x]. 
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Consider now the orbit of (x,Yq) under the action of the group G: 

(x, Y 0 )<^(xY 1 ,Y 0 )<^(xY 1 ,x 2 Y 0 )<-^(x, x 2 Y 0 )+^->(x, z 2 Yi)e-5->- 

(xY 0 , x^^ixYo, Y 1 )^(x, Fi). 

By construction, each pair (x',y r ) in this orbit cancels the kernel K(x,y). If more¬ 
over the series Q(x',y') is well-defined, then we derive from (8) that 

R(x') + S(y') = x'y', 

where we use the notation 

R(x) = t(Q(x,0) - Q(0,0)), S(y) =t(l + y)Q(0,y). (10) 

Since Yq is a power series in t, the pairs (x,Fo) and (x,x 2 lo) can obviously be 
substituted for (x,y) in Q{x,y). Given that xYq is symmetric in x and x, these 
pairs are derived from one another by replacing x by x. One has to be more careful 
with pairs involving Fi, since this series contains a term x/t. With the step set 
that we consider, and the non-negativity conditions imposed by the quadrant, one 
easily checks that each monomial x t yH n occurring in the series Q(x,y) satisfies 
n + i — j > n/2. Given that Yq = 0(£) and Fi = 0(l/f), this implies that the 
pair (xFo, Fi) and its companion (xFo, x 2 Fi) (obtained by replacing x by x) can be 
substituted for (x,y) in Q(x,y), so that Q(xYo,Yi) is a series in t with coefficients 
in Q[x, x]. We will not use the other pairs of the orbit, and the reader can check 
that they do not give well-defined series Q(x',y'). 

We thus obtain a total of four equations: 


R{x) + 5(F 0 ) = xY 0 

(11) 

R(xY 0 ) + S(Y 1 ) = x, 

(12) 

R[x) + 5(x 2 F 0 ) = xF 0 , 

(13) 

R{xYq) + 5(x 2 Fi) = x. 

(14) 


3.2. An equation relating R(x) and R(x ) 

We will now construct from the above system two identities that are symmetric 
in F 0 and Fi, and extract their non-negative parts in x. 

We first sum the first two equations, and subtract the last two: 

R(x) — S(x 2 Y 0 ) — S(x 2 Yi) + x = R(x) — S(Yq) — S'(Fi) + x. 

For two indeterminates u and v, the coefficient of t n in the series S(u) + S(v) is a 
symmetric polynomial in u and v. Applying Lemma 2 to this coefficient (for any n) 
shows that the right-hand side of the above identity is a series in t with coefficients 
in Q[x]. But the left-hand side is obtained by replacing x by x in the right-hand 
side. This implies that both sides are independent of x and equal to their constant 
term, that is, to —5(0) — 5(—1) (by (10) and Lemma 2 again). Finally, since S(y) 
is a multiple of (1 + y) (see (10)), this constant term is simply —5(0). 

We have thus obtained a new equation, 

5(F 0 ) + 5(Fi) = R(x) +x + 5(0). 

Combined with (11) and (9), it gives 

5(Fi) — xFi = f?(x) + R(x) + 2x — 1/t + x + 5(0). (15) 

For our second symmetric function of Yq and Fi, we take a product derived 
from (11) and (15): 

(5(Fq) — xlo)(5(Fi) — xFi) = — R(x) (R(x) + R(x) + 2x — 1/t + x + 5(0)). (16) 
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We want to extract the non-negative part in x. Let us focus first on the left-hand 
side. By Lemma 2, the term S(Yo)S(Yi) contributes 5(0)5(— 1), which is zero since 
S(y) is a multiple of (1 + y). Then x 2 YqY\ equals 1 by (9). We are left with the 
term — x(Lo5(Yi) + Yi5(Yo)). By Lemma 2, the factor between parentheses is a 
series in t with polynomial coefficients in x, and its constant term in x is — 5(0). 
But we also need to determine the coefficient of x in this series. Expanding S(y) 
in powers of y shows that 

YoSQT!) + Ei5(F 0 ) = (>o + yi)S(O) + Y 0 Y 1 F(Y 0 , Yi), 

for some series F(u, v) in t with symmetric coefficients in u and v. Since YqYi = x , 
the coefficient of x in YqS(Yi) + Y\ 5(Yq) is the same as in (Yq + Yi)5(0), namely 
5(0)/f. We can now extract the non-negative part from (16), and this allows us to 
express the non-negative part of P(x)P(x): 

1 + (x — l/t)5(0) = —P(x) 2 — [x-]P(x)P(x) — (2x — 1/t + x + S(0))R(x). 

(Recall that R(x) is a multiple of x.) Extracting the constant term in x gives 

1 - S(0)/t = ~[x°}R{x)R{x) - 2f?'(0). 

Since R(x)R{x) is symmetric in x and x, we can now reconstruct it: 

R(x)R(x) = [x-]R{x)R(x) + [x~]l?(a:)l?(ir) — [x°]i?(x)i?(x) 

= — R(x) 2 — (2x — 1/t + x + S'(0))i?(x) — R{x) 2 

- (2x - 1/t + x + S{0))R(x) - 1 - (x + x - l/t)S{0) + 2R'(0). 

That is, 

R(x) 2 + R(x)R(x) + R(x) 2 + (2x — 1/t + x + S(0))R(x) 

+ (2x - 1/t + x + S(0))R(x) = 2R\0) -(x + x- 1 /t)S(0) - 1. (17) 

3.3. An equation for R{x) only 

We cannot extract the positive part of the above equation explicitly because of 
the “hybrid” term R(x)R{x). However, the form of the first three terms suggests a 
multiplication by R(x) — R(x) to decouple the series in x from the series in x. More 
precisely, if we multiply (17) by R{x) — R(x) + x — x, we find a decoupled equation 

P(x) = P(x), (18) 

with 

P(x) = R(x) 3 + (5(0) + 3x — 1 /t)R(x) 2 

+ (2x 2 — x/t + x/t — x 2 — 2R'(0 ) + (2x — l/t)S(0)) R(x) 

- x 2 S( 0) + x (2i?'(0) + 5(0)/£ - 1). (19) 

Given that R(x) is a multiple of x, all terms in the expansion of P(x) have an 
exponent (of x) at least —1. But then (18) implies that P(x) is a symmetric 
Laurent polynomial in x, of degree 1 and valuation —1. More precisely, 

P(x) = [x°]P(x) + (x + x)[x]P(x) = [x°]P(x) + (x + x)[x]P(x). 

Thus, by expanding the expression (19) of P(x) in x at order 0, we find: 

P(x) = 2(x + x)P'(O) + P'(0)(25(0) - 1/t) + R"( 0). 
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Returning to (19), this gives 

R(x ) 3 + (5(0) + 3x — l/t)R(x ) 2 

+ (2x 2 — x/t + x/t — x 2 — 2R'(0) + (2x — l/f)5(0)) R{x) 

= R"{ 0) + R'(0)(25(0) + 2x-l/t) + xS(0)(x - l/t) + x. (20) 
Thus R(x) satisfies a cubic equation over Q(t, x , 5(0), i?'(0), 1?"(0)), 

Pol(R(x), 5(0), R'(0), R"(0),t, x) = 0, 

with 

Pol(xo, X 1; X 2 , X3, t, x) = Xq + (xi + 3x — 1 /t)XQ 

+ (2x 2 — x/t + x/t — x 2 — 2 x 2 + ( 2 x — l/t)x\) xq 

— X 3 — X 2 ( 2 x\ + 2 x — l/t) — xx\(x — l/t) — x. (21) 

In particular, if we prove that the one-variable series 5(0), R'( 0) and 1?"(0) are 
algebraic over Q(t), then (20) shows that R(x) is algebraic over Q(t,x). Observe, 
as an encouraging sign, that R(x) = t(Q(x, 0) — Q(0,0)) has degree (at most) 3 
over Q(t,x, 5(0), i?'(0), i?"(0)), in accordance with Theorem 1. Equation (20) is, 
in disguise, an ordinary discrete differential equation of degree 3 and order 3, as 
discussed in Section 2. 

3.4. The generalized quadratic method 

We have described in [15] how to study equations of the form 

Pol(R(x),Ti,..., A k ,t,x) = 0, (22) 

where Pol(xo, Xi,..., x*,, t, x) is a polynomial with (say) rational coefficients, R(x) = 
R(x; t) is a formal power series in t with coefficients in Q[x], and A±, ..., Ak are k 
auxiliary series depending on t only (in the above example, Pol is a Laurent poly¬ 
nomial in t and x, but this makes no difference). The strategy of [15] instructs us 
to look for power series X(t) = X satisfying 

d Pol 

— (R(X),A 1 ,...,A k ,t,X) = 0. (23) 

Indeed, by differentiating (22) with respect to x, we see that any such series also 
satisfies 

d Pol 

— (R(X),A 1 ,...,A k ,t,X) = 0, (24) 

and we thus obtain three polynomial equations, namely Eq. (22) written for x = X , 
Eqs. (23) and (24), that relate the (fc + 2) unknown series R(X), Ai,... ,A^ and X. 
If we can prove the existence of k distinct series Xi ,..., X satisfying (23), we will 
have 3 k equations between the 3k unknown series R(X 1 ),..., R(Xk), Ai ,..., Ak, 
Xi ,..., Xfc. If there is no redundancy in this system, we will have proved that each 
of the 3k unknown series is algebraic over Q (t). 

We apply this strategy to (20), with A\ = 5(0), A 2 = i?'(0) and A 3 = R"( 0). 
Equation (23) reads: 

3 R(X ) 2 + 2(5(0) + 3/X - 1 /t)R(X) 

+ 2/X 2 - l/(tX) + X/t -X 2 - 2R'(0) + (2/X - l/t)S{0) = 0. (25) 

Recall that R(x) := t(Q(x, 0) — Q(0,0)) and 5(0) := tQ{ 0,0) are multiples of t. 
Hence, once multiplied by tX 2 , this equation has the following form: 

X(1 -X)(l + X)=t Poll (Q{X, 0), R'(0), Q(0,0), t, X). 
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This shows that there exists exactly three series in t, denoted Xq, Xi and X 2 , that 
cancel (25). Their constant terms are respectively 0, 1 and —1. Due to the special 
form of our polynomial Pol (given by (21)), it is in fact simple to determine these 
series. Indeed, we observe that 
/ f) Pol f) Pol \ 

tx 2 ( — -h x 2 ——— ) = (1 — x)(l + x){ 2 tx 2 + 2 t — x)(xox + x\x + 1), (26) 

y OX 0 OX J 

and since the three series X,; cancel both partial derivatives of Pol, we conclude 
that X\ = 1, X'i = — 1, while Xo is the only formal power series in t that cancels 
(2 tX 2 +2 1 — X) (a Catalan-ish series, the exact expression of which we will not 
use). The fourth factor in (26) cannot vanish for xo = R(X) and x\ = 5(0) because 
of the factor t occurring in these series. 

Let Disc(x) be the discriminant of Pol(xo> 5(0), R'( 0), x) with respect to 

xq . It is obtained by eliminating xq between this polynomial and its derivative with 
respect to xq. Then Disc(Xi) = 0 for i = 0,1,2, and this gives three polynomial 
equations relating 5(0), i?'(0), R"{ 0), t and X 0 . In fact, we observe that Disc(x) 
is symmetric in x and x (this comes from the construction of (20)) and is thus 
a polynomial in s = x + x. So we can alternatively write that this polynomial 
vanishes at s = ±2 and s = 1/(2£). Systematically eliminating two of the three 
series 5(0), R'( 0) and R"(Q) in the resulting system of three equations yields an 
algebraic equation for each of them: for R'( 0) this equation has degree 4, and degree 
8 for the other two series. The equations for R'( 0) and 5(0) are shown below. From 
this, one can derive an algebraic equation for R(x) (using (20)), that is, for Q(x, 0), 
and then one for S(y) = t( 1 + y)Q(0,y), as explained below. Finally, the original 
functional equation (8) proves the algebraicity of Q(x,y). 

3.5. Rational parametrization 


To avoid handling big polynomials, it is convenient to use rational parametriza- 
tions of all these series. To begin with, the equation satisfied by I?'(0) is: 


729t 6 i?'(0) 4 + 243f 4 (4t 2 + l)R'(O) 3 - 27t 2 (14t 4 + 19f 2 - l)i?'(0) 2 

- (20t 2 - l)(7t 2 - 6t + 1)(7 1 2 + 6 1 + l)R'{0) - t 2 (343f 4 - 37 1 2 + 1) = 0. 


It has genus 0 —so Maple tells us— and can be parametrized by introducing the 
unique formal power series T = T(t) in t with constant term 1 satisfying (3). Then 


R\ 0 ) 


(T — 1) (21 - 6T + T 2 ) 
(T + 3) 3 


Rational parametrizations can be computed via the Maple command parametrization. 
The one that we originally obtained has a slightly different form, see (4). 

The equation satisfied by 5(0) is: 


27f 7 5(0) 8 + 108f 6 5(0) 7 + 189f 5 5(0) 6 + 189t 4 5(0) 5 - 9f 3 (32f 4 + 28t 2 - 13)5(0) 4 
- 9f 2 (64t 4 + 56f 2 - 5)5(0) 3 - 2t(256f 6 - 312f 4 + 156f 2 - 5)5(0) 2 
- (321 2 - l)(4t 2 -6 1+ 1)(4 1 2 + 6 t + 1)5(0) - t(256t 6 + 576t 4 - 48t 2 + 1) = 0. 


Since 5(0) = tQ( 0, 0), this gives an equation for Q(0,0), involving only even powers 
of t. If we replace t 2 by its rational expression in T derived from (3), this equation of 
degree 8 in Q(0,0) factors into a quadratic term and one of degree 6. Injecting the 
first few coefficients of <5(0, 0) shows that the term that vanishes is the quadratic 
one. So Q( 0,0) has degree 2 over Q(T), and can be written in terms of Z = y/T as 
stated in Theorem 1: 


Q(o,o) 


32 Z 3 (3 + 3Z -3Z 2 + Z 3 ) 
(l + Z)(Z 2 +3 ) 3 
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The series R"(0)/t is also found to have a rational expression in Z: 

R-( 0) = MX) = 1 024 t ^ Z - 1 )< 1 ; i 2 + Z + ^^I- 2Z5 + Z6 ) . 

We now go back to the equation (20) satisfied by R(x ) = t(Q(x, 0) — <9(0,0)). 
Injecting the above expressions of S'(O) = tQ( 0,0), I?'(0) and R"( 0) in it gives 
a cubic equation for <9(:r,0) over Q (t,Z,x). We now consider Q(xt, 0) instead of 
Q(x, 0), since it is an even series in t. Using (3) to express t 2 in terms of T = Z 2 , 
we see that this series is cubic over Q(Z,x). 

We finally construct an equation for Q(0,y ) thanks to the kernel equation (11). 
Written for xt instead of a;, it reads 


t(Q(xt, 0) - Q(0, 0)) + t( 1 + y)Q( 0, y) = xty, 

with y = Yo(xt). Thus the equation we have obtained for Q(xt,0) gives a cubic 
equation for <9(0, y), and eliminating x between this equation and K(xt,y) = 0 
gives a (still cubic) equation for <9(0, y) over Q (Z, y), valid for a generic value of y. 

The equations for Q(xt,0) and <9(0, y) take a few Maple lines each (four or 
five), and we do not give them explicitly. They can be seen in the accompanying 
Maple session. They factor if we express x and y in terms of van Hoeij’s series U 
and V, as described in Theorem 1. Our alternative parametrization, given below 
the theorem, was obtained thanks to the parametrization command in Maple. 
One can use this command with Z as a parameter, but it is faster to compute a 
parametrization for a few values of Z (e.g. Z = 10, Z = 100) and then reconstruct 
a generic one. □ 


4. More algebraic models 

It is natural to ask when the symmetric functions of the roots of the kernel 
are polynomials in x, since this property plays a crucial role in our proof. Let us 
consider each of the 23 models with small steps and a finite group [16], and denote 
by Yq and Y\ the roots of the kernel 

R{x,y) = i -t ^2 xl y°- 

It is not hard to see that lo + Y\ and YqY\ are polynomials in x if and only if the 
only step of the form (j, 1) is (1,1). We find 4 models having this property. 

• The first one is Gessel’s model S = {—►,/*,•<—,,/}, which we have just 
solved. 

• Then comes Kreweras’ model S = {■<—,4-,/’}, which is also algebraic, and 
was solved in [11] using the same principles as in this paper (see [12] for a 
variant). The solution is simpler than in Gessel’s case in two aspects: first, 
the diagonal symmetry implies that Q(x, 0) = <9(0, x), so that we have only 
one unknown series R(x), not R(x) and S(y) as before; then, and more 
importantly, the equation obtained by forming a symmetric function of To 
and Y\ only involves R(x ), not R(x). 

• Finally, we also have the models {•<—, \, /*} and {•<—, -A, \, /'}, which are 
known to be D-finite [16] but transcendental (this can be derived from [10, 
Thm. 4]). 

More recently, Kauers and Yatchak initiated a study of walks in the quadrant 
with multiple steps [32]. Such walks naturally arise, after projection, in the study 
of 3D walks confined to the first octant [6]. In particular, it was proved in [6], 
using computer algebra, that the model —>, —>, /*} is algebraic (note the 

double East step). This was generalized by Kauers and Yatchak [32], who proved 
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that algebraicity persists if one includes a South step, with arbitrary multiplicity A. 
For this model, the symmetric functions of the roots of the kernel are polynomials 
in x, and we can solve it using the tools of this paper. In fact, the proof is only 
marginally more difficult than in Kreweras’ case: there are two unknown series 
R(x) and S(y), but the equations obtained by forming symmetric functions of Yq 
and Yi only involve R(x), not R(x). 

Let us briefly sketch the main steps of the solution. The basic equation reads: 


K(x,y)xyQ(x,y) = xy - R(x) - S(y), (27) 

where K(x , y) = 1 — t(x + xy + Ay + xy + 2x + xy) is the kernel, 

R(x) = t(l + Xx + x 2 )Q{x, 0) — tQ(0, 0) and S(y) = t(l + y)Q(0,y). (28) 


The group of this walk has order 6, and the orbit of (x, To) contains exactly 4 pairs 
(x',y') for which Q(x',y') is well defined: 


{x,Yq), 


\ 1 T At 



( t(l + Y 0 ) 
\ 1 + At 



( ^(f + Tp) 

\ 1 + At 


c/t + Xx — 1 ) . 


Each such pair gives an equation x'y' = R(x') + S(y r ). The first equation expresses 
S(Yq) in terms of R(x), and a combination of the last two expresses S(Yi) in terms 
of S{x/t + \x — 1). By extracting the positive part in the resulting expression of 
S(Yq) + S'(Yi), we conclude that 


S'(Yo) + 5 , (Ti) = x and x = S(x/t + Xx — 1) — R(x). (29) 


From this we form a second symmetric function of To and Y\: 

(5(T 0 ) - xY 0 )(S(Y 1 ) - xYi) = -R(x) {R{x) +2x + 2x- 1/t). 


Extracting the positive part in x gives 

Xx + x 2 = —R(x) (R(x) + 2x + 2x — 1/t) + 2R'(0). 


This discrete differential equation can be solved using the method of Section 3.4. 
One obtains a cubic equation for i?'(0), and one of degree 6 for R(x). Finally, 
saying that R(x) equals —tQ{ 0,0) modulo (1 + A.t + x 2 ) (see (28)) gives, by taking 
the above equation modulo (1 + Aa; + x 2 ), 

t 2 Q{ 0,0) 2 + (2A t + 1)Q(0,0) = 2i?'(0) + 1, 

and proves the algebraicity of Q(0,0). One can then compute an equation for 
S(y) using the second part of (29). The algebraicity of Q(x,y) follows using the 
equation (27) we started from. □ 


To complete our guided tour of (conjecturally) algebraic models, let us mention 
three models that were discovered by Kauers and Yatchak [32]: 

2 12 



The weights indicate multiplicities. All three models have a group of order 10. We 
can solve the first two using half-orbit sums, as in [16, Sec. 6], and we find them to 
be algebraic indeed. Regarding the third model, its algebraicity should follow from 
a complex analysis approach [44], as for Gessel’s model in [8]. 

Finally, let us recall that another natural way of weighting steps comes from the 
study of Markov chains in the quadrant [23] . In this setting each step is weighted by 
its probability. It may be possible to construct a Markov chain with Gessel’s steps 
and stationary distribution given by an algebraic generating function, provable by 
the tools of this paper. This would parallel the case of Kreweras’ walks: they first 
appeared as a counting problem in 1965 [35], and then 20 years later, independently, 
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as a Markov chain with algebraic stationary distribution [26]. Many years later, the 
algebraicity of both problems was proved by a common elementary approach [12], 
using ideas that we have extended in this paper. 
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